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ABSTRACT 


The thermal signature of a surface vehicle, ‘“SPRUANCE"” 
type destroyer, quantified by the signal to noise ratio of a 
thermal seeking missile head is parametrically investigated. 
The effects of the ship internal temperature, ship body 
paint emissivity, sky condition, sun elevation angle above 
the horizon, atmospheric profile, missile optics and flight 
altitude are examined in detail. Results show that both the 
ship body temperature and signal to noise ratio increase as 
the incident solar energy and the ship body paint emissivity 
increase and that the signal to noise ratio appears a peak 
for sun elevation angle in the range of 40° to 60°. Moreover 


the signal to noise ratio increases as the missile flight 
altitude decreases and keeping the other parameters constant 
higher values are found for the Midlatitude Summer atmos- 
pheric profile than for the Tropical atmosphere. 
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I. INTRODUCTION 


The basic function of a thermal-seeking-head missile is 
to track on a surface vehicle using the vehicle's radiant 
intensity. This operation 1s based on the fact that a moving 
vehicle shows a generally higher surface temperature that 
the environment. Reasons for this are the existence of 
operating engines, resulting 1n both increased engine room 


temperature and “hot” combustion products exiting into the 
atmosphere from the stack, the operation of a number of 
devices which result in heat generation, the requirement for 
certain compartment temperatures comfortable to the crew, 
and the incident solar energy on the painted metallic 
vehicle surfaces. 

Exact information about the heat seeking missile charac- 
teristics is not available because of their classification. 
Hence our study will be a parametric one in order to examine 
the relative influence of a number of parameters on the 
vehicle thermal signature which is quantified by the signal 
to noise ratio of the thermal detecor head. The range of the 
parameter values will be chosen using both experience and 
unclassified information. 

Exceptional cases such as operation of the vehicle close 
to the coast or other vehicles, which could have similar or 
higher temperatures, will not be considered. At this point 
we must realize that the thermal interaction of the vehicle 
with the environment is a very complicated phenomenon, and 
its exact determination for the different parts of the ship 
body is outside of the scope of the present work. Instead we 
will consider a simplified case of uniform interaction for 
all the ship parts with the environment (convection, radia- 


tion, solar absorption energy) and will limit the source for 
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the thermal head detector enviromental noise only to that of 
the sea water. 

The main procedure steps of the present study will] be 

To create a simplified model for the calculation of the 
vehicle body temperature and to show what parameters affect 
it and in what direction. 

To determine the “hottest” vehicle area and combining 
with the operating characteristics of the missile to calcu- 
late the dimensions of the tracked area of the vehicle. 

To calculate the total radiant intensity of the tracked 
area of the vehicle in a certain wavelength range, where the 
thermal detector head is operating. 

To calculate’ the atmospheric transmittance, using 
LOWTRAN 6, for each given combination of the parameters of 
our interest in the given detector wavelength range and 
consequently the radiant heat flux detected by the thermal 
head. 

To create an appropriate detector model that could show 
optimum operation for the conditions of our problem, and 
finally to calculate the creating signal-to-noise ratio 
which quantifies the vehicle thermal signature and describes 
the detector operational quality. 

To discuss the relative influence of all the pertinent 
parameters and to propose certain considerations about both 
the vehicle and the missile. 

Without expecting significant quantitative importance 
for our results, due to the simplifications that have been 
made, their qualitative importance must be more remarkable. 

Parameters that will play significant roles in our study 
will be 

vehicle internal temperature, 

vehicle body paint emissivity, 

sky condition, 

sun elevation angle above the horizon, 
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atmospheric profile, 
missile flight altitude, 

instantaneous Field of View (I.F.V.) created by the 
detector'’s optic system, 

detector material, 

The present work consists of five Chapters. In the first 
Chapter we state the assumptions of the problem, we calcu- 
late the vehicle body temperature, and we determine its 
“hottest” area which is tracked by the detector. In the 
second Chapter we use the LOWTRAN 6 program in order to 
determine the atmospheric transmittance and in the third 
Chapter we calculate the total vehicle and background 
radiant intensity incident on the tracker aperture. In the 
fourth Chapter we calculate the detector signal-to-noise 
ratio for all the possible combinations of the parameters, 
and finally in the fifth Chapter we discuss the relative 
importance of each parameter and we Suggest ways that could 
improve both the vehicle thermal signature and the detector 


operation. 


Il. SHIP MODELING 


In order to create an exact ship model, it is necessary 
to acquire a significant amount of empirical data. For the 
present study, and for several reasons, the available data 
are limited, and are based only on the unclassified litera- 
ture given inthe first six references, with relatively 
reasonable assumptions. 


A. ASSUMPTIONS 


All compartments of the ship have a uniform temperature 
of 20°C=293.15K except the main engine rooms 1 and 2, which 
have a temperature of 30°C=303.15K. 

The stack exit temperature is uniform and we study the 
heat radiation from the exit plane neglecting the plume 
radiation. Reasons for doing this are the non-uniformity of 
the plume geometry and temperature distribution, and the 
fact that it is cooled significantly in a relatively short 
distance from the exit plane {[Ref. 3]. We choose as stack 
exit temperature 306°F=152.2°C=425.4K 

In order to account for the heat convection from the 
ship to the air we will use the relative wind-speed of 30 
Knots. For this case we will have forced convection with a 
heat transfer coefficient in the range of 5 to 10 (BTU/H 
FT?R) or equivalently 28.4 to 56.8 (W/M?K). We will not 
account for the heat leakage from the ship body to the sea 
water or for the heat loss due to the vaporization effect on 
the wet surfaces above sea level. For this reason and in 
order to be closer to reality we choose for our heat convec- 
tion losses calculation the value of 10 (BTU/H FT?R) or 56.8 
(W/M’?K) as representative for the convection heat transfer 


coefficient. 
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The present study will be in Mediterranean atmospheric 
conditions and more specifically for the area of the Aegean 
Sea. We will use two atmospheric profiles i.e. Midlatitude 
Summer for months October through April and Tropical for 
months May through September. The corresponding ambient 
temperatures are 21°C=294.2K and 26.6°C=299.7K respectively 
(Ref. 6]. 


B. GEOMETRY AND MATERIAL OF THE SHIP 


In Figures 2.1 and 2.2 you will find the ship side and 
top view respectively with the appropriate dimensions. 

The material of the ship lower body construction is 
constructional carbon steel with thickness 13.5mm = and 
thermal conductivity 54 (W/M K) (Ref. 9]. 

The material of ship body superstructure is duralaminum 
with thickness 3mm and thermal conductivity 164 (W/M K). 

These small values of the thickness and the high values 
of the thermal conductivities for both materials result in 
negligible temperature gradients inside the materials. Hence 
it is reasonable to assume a uniform temperature throughout 
the thickness of both materials. 

Both structures are painted with gray paint i.e. the 
emissivity and absorptivity will be spectrally independent 
and equal between them. A typical value for oil paint emis- 
sivity - absorptivity is 0.94 [Ref. 10:p.07] and for TiQ, 
gray paint is 0.87 [Ref. 4:p.278] both for temperatures 
300K. From these values we observe that the available data 


give us paint emissivity - absorptivity values in the range 


of 0.90. Based on the facts that these values play a signif- 


icant role in the determination of the ship body temperature 
and that intensive research resulted in the production of 
new paints with half the values of the already existing ones 
(Ref. 7] we will use two (2) different values i.e. 0.80, 


oO 


0.55. Therefore we cover the useful range of the existing 
and the newly developed paints. At this point we must note 
that we don't have any available information about selec- 
tively absorbing paints that could exist and give the 
minimum possible thermal signature and our assumption about 
the equality between emissivity and absorptivity is inevi- 
table as a first approximation. Moreover the solution of the 
problem for two (2) different emissivities constitutes a 
good research point for the relative influence of this 
factor. 

The emissivity of the interior surface of the ship body 
will be approximated as that of stainless steel which is 
0.44 (Ref. 4:p.278]. 


C. SOLAR HEAT FLUX IN THE SEA LEVEL SURFACE 


In Figure 2.3 you will find the illuminance levels on 
the surface of the earth due to the sun for three different 
sky conditions i.e. Unobscured Sun, Sun with Light Clouds 
and Sun with Heavy Storm Clouds as a function of the eleva- 
tion angle above the horizon. These three sky conditions and 
the elevation angle above the horizon will be two more 
parameters in the solution of our problem. 

Since the given values in Figure 2.3 are in photometric 
units (Lux) and we want to convert into useful radiometric 
units (W/M?) we have to determine a conversion factor. 

We know that the maximum incident solar flux on the 
earth surface is approximately 900 (W/M?) (Ref. 10:p.11]. 
This value coresponds to 1.15°10*% (Lux) and therefore the 
conversion factor will be: 


CF 


un 


900 / 1.15°10° = 5.826°10% (W/M?Lux ) 
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Using this conversion factor and the values of Figure 
2.3 we chtate the mormal Component of the solar heat flux in 
tw “4 ) oof 2 surface at sea level. Table 1 demonstrates the 
above Vblues. 

In order to calculate the values of the solar heat flux 
fallimy oblaicueisv ot as surface we must aiantroduce the inci- 
aetit Stese FP ardomeltsply the values of Table 1 by cos@. 

For a vertical tace the angle 0 1s identical with the 
esevatibon anepbe above the horizon. For a horizontal face 
Loe the te. F it our case, the angle @ 1s the complement of 


elevat ton ateee above the horicon. 


D. CALCECLATION GF SHIP BOY TEMPERATURE 
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Figure 2. Santrol Volume for Determination of Ship 


Body Temperature 


In Figure 2.4 you will find the appropriate control 
volume of the ship body and the pertinent energy flows that 


are related with ait. For simplicity the control volume has 
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unit area exposed to both the environment and the ship 
interior. More specifically: 
paint emissivity 
paint absorptivity (= a.) 
emissivity of the interior surface. (€m= 0.44 ) 
ambient temperature 
ship compartment temperature 
normal component given by Table 2 
net heat exchange of control volume with the 
environment 
net heat exchange of control volume with the 
interior of the ship 
h : convection heat transfer coefficient with the air 
The absorbed solar heat flux is Ey*a,q. The net heat 
exchange with the environment is due to heat radiation and 
heat convection i.e. 


Q,740( T*-Ty )+h( T-T.) (2.1) 


The net heat exchange with the interior is due only to heat 
radiation i.e. 


Qn=er( T*-Ty ) (2.2) 


Energy balance of the control volume under thermal equilib- 


rium gives: 


Ese? a, =60( T*-To )*h(T-T)*@o( T'-TS ) 


F(T)=(€,e€q)OT* *hT-o(eTy +Gals )-HTy Esq? ap =O 


a ils -s a 


Equation 2.4 must be solved numerically. We will use the 
Steffensen Recursive Method [Ref. 11]. The idea of this 
method is that starting with an x close enough to the exact 
solution S of the equation F(x)=0 then the sequence: 


Kee Xu LER (Xa)! EC xae £ (xn) )- F(xn)] (2.5) 


will converge very fast to the exact solution S. 

In our case the equation is F(T)=0 and we are looking 
for the value of T that satisfies it. 

In order to obtain an initial value so that our sequence 
will converge we use the value: 


Tin= [ (Egn® ae *8,0T,* * E6T in ) / (+m) o) 4 (2.6) 


The above initial value of T does not account for the 
heat convection and for our purpose gives acceptable initial 
values. Finally we stop the recursive formula when |TU-Ts | 
is less than or equal to 10*. 

In Appendix A, Tables5 through 18 demonstrates the ship 
body temperature for both horizontal and vertical faces, for 
all the combinations of two (2) atmospheric profiles, three 
(3) different sky conditions, two (2) internal (Ti) tempera- 
tures and two (2) paint emissivities. 


E. DETERMINATION OF THE SHIP “HOTTEST” AREA 


From the calculations in section D,we observe that the 
higher ship body temperatures happen to be in the region of 
the two engine rooms (higher internal temperature). 

From the ship geometry,it is obvious that the dimensions 
of the engine room No 1 are greater than those of engine 
room No 2. Therefore the total radiant intensity (W/Sr) of 
engine room No 1 will be greater since for both regions we 
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have the same normal components of radiance (W/cm?Sr). In 
addition to that to each engine room there corresponds the 
Same stack geometry also going a ‘hot’ region too. 

Under this simplified assumption we conclude that for a 
heat seeking missile head the area of engine room No 1 from 
the stack exit all the way down to sea level is the most 
significant region. | 

At this point it is necessary to couple the ship with 
some characteristics of the missile detector operation. 

Without getting involved in details of the detector 
operation, which will be studied in later sections, we will 
determine the detector Instantaneous Field of View (I.F.V.), 
so that we will know exactly the dimensions of the ‘hottest’ 
area which will be tracked by the detector and therefore we 
will be able to calculate the amount of radiant energy that 
will be transmitted to the missile. 


F. DETERMINATION OF THE DETECTOR INSTANTANEOUS FIELD OF 
VIEW 


The following range of data is based on our experience 
and research, aiming to calculate a reasonable and useful 
I.F.V. In Figure 2.5 you will find the corresponding 
geometry. 

The range of the missile flight altitude will be from 
50m to 100m above sea level. Since we want to reduce as much 
as possible the background noise, the detector must have an 
inclination with respect to the horizontal. We choose the 
range for this inclination (angle @) from 2° degrees to 5° 
degrees. 

From general cosideration for the detectors [Ref. 4] we 
know that I.F.V. is a function of the ‘hottest’ ship area 
that the detector has to track. 
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For the creation of I.F.V. it is necessary that the 


detector accepts radiation from inside an area which is 


related to an angle 26 generated by the detector mirror 
system as in Figure 2.5 . From the same figure we conclude 
the following: 
@ = 90°-¢ -& 
BC = AB + tan@ 
AC (AB? + BC?) 
CF AC*sin6& 
CG 2CF = 2AC*sinS = SIDE 
AF CG/2tan6 = PATH 
Where: 
CG represents the side of the I.F.V. 
AF represents the distance or path length between detector 
and the I.F.V. 
AB =H: flight altitude 
For simplicity we will consider that the shape of I.F.V. 
is square with side equal to CG of Figure 2.5 
In Appendix B Tables 19 through 29 demonstrate the 
values of PATH and SIDE for all the possible combinations of 
g¢, © and H. Program SHIP 1, in Appendix E, has been used for 
the calculation of the above values. 
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G. DETERMINATION OF THE DIMENSIONS OF THE TRACKED SHIP AREA 


In the previous section we have shown that the detector 


will track on the area of engine room No l. 
In Figure 2.6 you will find a three-dimensional sketch 


for the engine room No 1 area. Before we start the calcula- 
tions of the different surfaces we will explain why the side 
‘ of I.F.V. is greater than the width of the engine room No 1 
which is Wi=},620cm. We consider as useful that the I.F.V. 


oor will include the whole area between sea level and the stack 


exit plane. The total height of this dimension is: 
= Hit Hyg+ Hat Hg =720+900+360+180=2,160cm 

From the Tables 13 to 24, we choose for each flight 
constant 6=0.56° SIDE that is 


the minimum 


altitude and for 


greater than or equal to 2.160cm 
Table 2 demonstrates the results of this choice 


ney 


Based on Figure 2.6 we will now determine the character- 


istics of each individual surface. 
A,= Wie Hy 
s Vertical face | 
oa Tw= 303.15K 


Au= (SIDE - W, )°H, 
Vertical face 
Tw= 293.15K 


Ag= Wie Ha 
Horizontal face 
Tin= 303.15K 


Aas (SIDE ng W, )°H, 
Horizontal face 
T= 293.15K 


gS arian gig aoe 
Se Pee Gt oe 


Beary ,,3S99390H,, aya JO YIAAS [euOTSuUaUWTG-saIY] 9°Z7 aaindty 
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ee TABLE 2 
ot Final Choice for the Dimensions of I.F.V. 
Te 
HEIGHT SIDE PATH DELTA PHI 
FLIGHT (m) m) 
“$0 (m) | 21.88 | 1119.34 | 0.56° | 2.0° 
60 (m) | 21.97 | 1123.89 | 0.56° | 2.5° 
. “70 (m) | 22.04 | 1127.25 | 0.56" | | 3.0° 
“80 (m) | 22.09 | 1129.84 | 0.56° | 3.5° 
= "90 (m) | 22.13 | 1131.96 | 0.56° | 4.0° 
a “100 (m) | 22.16 | 1133.72 | 0.56° | 4.5e 
x OE se VE se eg Se Ae oa ee eS oe os Sa ee ee he ee 
re As= SIDE ¢ Hs 
Vertical face 
Ta= 293.15K 
7 
‘ Aq= Wy? Hy 
Vertical face : 
Tw = 303.15K 
At this point we must notice that surface As corresponds 
oS) to the stack and,since we don't have information about the 
| insulation of this structure, we approximate its temperature 
with that of the surface with T = 303.15°K 
; Ag= 4.43°10* cm? 
9. Horizontal face 
- 2 Tn= 425.4K 
‘ Ay= SIDE*H,- (1/4) °W2 


Horizontal face 
Tw= 293 .LSK 


a For the order of magnitude of PATH i.e. the distance 

- between missile and ship we will approximate the thermal 
emission of each individual surface as point source emis- 

ok sion. These "point sources" will be: 

aoe (Au+ Au), (Aa+ Ana), As, Aus Ags Acs Az, 


- with corresponding mean altitudes from sea level: 
Hew= Hew® Hy /2 

= Hay? Hac Ha 

oe Hay? Hy + Hy /2 


Ha= H, + Hy + Hy /2 
Hae Hy + Hy*+ Hat Hy /2 
Ha= H; ¢ H, + Hy, + He 


Hae H,+ Hs 
The data for the above calculations are: 
W. = 1,620cm 
H,= 720cm 
H,= 210cm 
H,= 900cm ; 
H,= 1,260cm 
oe W, = 900cm ’ 
oe H,= 360cm 


: Hg = 180cm 


III. ATMOSPHERIC TRANSMITTANCE USING LOWTRAN 6 


A. DESCRIPTION OF LOWTRAN 6 


code that has 


LOWTRAN 6 is the latest issue of LOWTRAN 
been designed to calculate the transmittance of the atmos- 


phere over spectral bands from 0.25pm to 28.5pm along paths 


of various inclinations and, of course, through different 


atmospheric profiles and parameters. 
In order to use LOWTRAN 6 program we have to fill four 


(4) cards (groups) of input data. 


Cards 1 and 2 require the significant meteorological 


parameters for the atmospheric profile. Card 3 requires the 


Ld geometry of the problem i.e. initial and final altitude and 


the path length. Card 4 is related to the wavelength range 


and limits. 


B. DATA USED 


For card 1 we used data for the Mediterranean Sea and 


more specifically for the Aegean Sea area and the region 


r) from Crete to North Africa. 
For this area we have two different atmospheric 
profiles, namely Midlatitude Summer for months October 


through April, and Tropical for months May through September 
[Ref. 6]. 
For card 2 we used again data from [Ref. 6] for each 


month individually and for card 3 we used all the possible 
final altitude and PATH from Table 2 
altitudes of each 


combinations between 
for the 
ie “point source" of the "hottest" ship area as given in 
: Section G of Chapter II. 


and initial altitude the mean 


a i 
: ml 


At this point we have to choose one more parameter that 
will be needed for our problem : the wavelength range of 
the detector operation. We decide that this range will be, 
in first approximation, in the Sum window (i.e. AMW=4.20um - 
wz5.18pm) [Ref. 4:ch 5, p.89]. 

Lowtran 6 has the advantange of giving us, besides the 
average transmittance between two limits, the spectral 
transmittance as well. From this spectral distribution we 
can reject certain values of very low spectral transmittance 
and therefore the wavelength range can be narrower. This 


peculiarity will have an influence on card 4. 


C. RESULTS 


Running LOWTRAN 6 for all the flight altitudes of the 
missile and for both Midlatitude and Tropical atmospheric 
profiles, we observed that, for wavelengths less than 
4.435ym and greater than 5.0pm, the spectral transmittance 
generaly shows values less than about 0.10. Deciding to 
reject these low values we choose as the final spectral 
range of the detector operation to be from 4.435um_ to 
5.0Oum 

In Figure 3.1 you will find the spectral transmittance 
in the above range for both atmospheric profiles. 


Using the final wavelength range we calculated the 


average transmittance for all the possible combinations 
between atmospheric profile, missile flight altitude and 
individual "point sources" of the ship's tracked area. 
Tables 23 to 46 in Appendix C demonstrate both the data and 
total average transmittance for each month. 

From the data of the above tables, it is obvious that 
the influences of both missile flight altitude (from 50m to 
100m) and the relative differences between the individual 


“point source" altitudes (4m to 16m) are very weak on the 


average transmittance. Averaging in both directions in this 
two dimensional table results in a unique representative 
value for each month. Checking these average values we 
observe again a relatively small difference between months 
belonging in the same atmospheric profile. Since a unique 
value of average transmittance of each atmospheric profile 
would be very convenient for the calculations to follow, we 
average over several months yielding the following final 
values: 

Midlatitude summer : 0.4711 

Tropical: 0.4259 
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IV. CALCULATION OF SHIP AND BACKGROUND RADIANCE 


A. NORMAL COMPONENT OF SHIP RADIANCE DUE TO THERMAL 
RADIATION 


The purpose of this calculation is to find, by inte- 
grating the Plank’s radiation law, the radiance in the given 
wavelength range. 


From Plank’s radiation law we know that: 


M,(T) ={2mcth a Je[1 (em -1)}  (Wocm? cm) 
or 

M (T) =Ga‘ (es 2k) 
where: 

c, 2 2we*h=3.7415¢L0%* (Wecm? ) 

ci, zhe k-1.4388 (cmek) 
The above formula is in the form of: 

M (T)=Cox™(e™ -1) 

Integrating this between the band limits y,, dy, we will 

find the in-band energy radiance (Ref. 4:p. 19]. 


Therefore: 


»3 
Ly: fem fMy(T)dy= (Cie CyjeTts2, (W/cm? *Sr) 
Where: me 
E21 S ems [ (mex)? 3 (mex)? +6 (mex)? + 6mex-6] eee 
meh a 
and x=C,/,AT 
E-Swe [ (mex )'+3(mex )*+6mex-6] - 
etin* ( (mex )'+¥imex ) + 6mex-6] 
In the above relationsythe temperature T is taken from 


Tables 5, 18 im Appendix A. 
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B. NORMAL COMPONENT OF SHIP RADIANCE DUE TO REFLECTED SOLAR 
ENERGY 


We have already calculated the normal component of solar 
irradiance E incident on each ship surface. In the energy 
balance calculation we considered that the solar’ energy 
absorbed by the surface is © ae. Obviously the amount 
(l-a.)*E,is reflected. 

Defining as D(0-A) the percentage of the solar constant 
associated with wavelengths shorter than A we take 
[Ref. 4:p. 17]. 

D(0-4.435)=99.3075454% 
D(0-5.000)=99.511500% 
Therefore 

AD=D(0-5.000)-D(0-4.435 )=0.2039546% 

In conclusion we have that the normal reflected solar 
energy per m* from 41=4.435pm to %=5.00pm is: 

(1-ap )*AD* Eye (1L-a, )°2.04°10° °E,, (W/m?) 

In order to determine the normal component of radiance 
due to reflection we must divide the above value by tT and in 
order to transform into cm? we must divide by 10°. Finally 
the normal component of radiance due to reflection in 
(W/cm?*Sr) from =4.435pm to %=5.00pm will be: 

Laaa= (1-ae )°2.04°10° *E,/m e1l0*=(l-a, )*4.08°10%*E,, (W/cm?eSr) 


C. NORMAL COMPONENT OF STACK EXIT PLANE RADIANCE 


So far we know that the stack exit plane area is 
4.43°10* cm? and that the corresponding uniform temperature 
is 425.4K. 

For simplicity we consider that the greater part of the 
exaust gas is CO,. Since CO, is a selective emitter we need 
to carry out the approximation shown in Figure 4.1 which 
partially accounts for atmospheric absorption and approxi- 
mates the peak with a constant emissivity over a limited 
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Figure 4.1 The 4.4um Emission Band of COqg[Ref. 12] 
range. Taking the band to be from %=4.35um to ”2=4.55um 
© and the emissivity 0.5 we find for the radiance: 


L eof tieda =(goT* /m VfMada/ [Ms daz (oT! /)* (0.05) 
=(0.5)°(5.67°10™ )°(425.4)"%*(0.05)/7 
=1.478°10% (W/cm?*Sr) 
¢ Since our band of interest overlaps with half of the 
4.4m emission band of CO, the above value must be divided 
by 2. Therefore: 
Pees L=7.39°10*% (W/cm?*Sr) 
* 
D. TOTAL SHIP RADIANT INTENSITY IN THE DIRECTION FROM SHIP 
TO MISSILE 


As mentioned in Section G of Chapter II,we will consider 
the surfaces (A,+Ay), (AatAa), Ag, Ag, Ag, Ag, Ag, aS “point 
sources". For each source we know the corresponding radiance 
due to thermal radiation and reflection i.e. 

Lay =La. +Las, 
Lr =Laan* Laaan 
Lara =Lay + Laaa 
Lraa = Laat Lagan 


Las =Lasy+Laaas 
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= Lon, +Learn, 

= Lays + Laaas 

=Lay +tLeass 
For Lre we have the unique value given in section C of 
chapter III 


In order to obtain the normal radiant intensity (W/Sr) 


of each surface we must multiply the corresponding values of 
each area and radiance. Doing this we take: 
Tye =Lm eA, +L? Au 
ZLreAy +Lan® Aa 
=Las®Ay 
=Lru* Ay 
=Las*As 
=Lret Ae 
Iyn =Li*Ay 
The important magnitude for us is the radiant intensity 


in the direction of source - missile. 


a tlbbateds edatdddtiad 


Figure 4.2 Relative Position Between Surface- Missile 


In Figure 4.2 you will find a sketch representing the 
relative position between the surface and the missile. Based 
on that we conclude this the radiant intensity will be: 


I=Iynecosx for a vertical face 
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I=In*sinx for a horizontal face 
where x=sin! (H-Hs»/ PATH) 


The resulting relations will be: 


I, =Iw*cosx 
Ia =Iwesinx 
Is =Iyn* cosx 
Ig =Iyy* cosx 
Ig =Isn® COSX 
I, =Laesinx 


I, =e sinx 


E. BACKGROUND RADIANT INTENSITY 


In the ideal case in which the detector tracks perfectly 
on the ships "hottest" area and the square I.F.V. will fit 
perfectly on that the detector will still see some areas 
outside of the ship. This is due to the fact that the upper 
construction of the ship has not a perfect rectangular 
shape. Therefore the detector will receive signal from an 
area of the sea that radiates behind the ship. Since this 
Signal constitutes noise for the detector, it has important 
meaning. 

In Figure 4.3 you will find the geometry from which we 
will determine the background radiant intensity. We must 
notice that the sketch in Figure 4.3 has been magnified for 
convenience in the ship region. In other words angle (98+6) 
is very close to 90° and the side of I.F.V. AE is ASE is 
practically equal to the ship dimension ABC. 

AB represents the height of surface A; [Figure 2.6] from 
the sea level or AB=H,+H, 

BC then represents the path through which the detector 
will see the sea area FG past the stack. 
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e The width of this area is (SIDE-W,) which represents the 
& free horizontal dimension over the surface Aj. 
r. , From the geometry of Figure 4.3 we have: 

. DS=AS-SIDE 

be @, =tan* (DS/PATH) 

ao Q, =6-O. 

co MA=PATH/cos8 

ea ; @=cos (H/MA) 

an PF=Hetan(@+6+, ) 

= PG=H*tan(@+26) 

7 FG=PG-PF 

O, =O+6+Q, +02 /2 

MK=HetanQ, 

ae ; AREA=FG* [SIDE-W, ] *cos@, 

CG Table 3 demonstrates the result of the above calcula- 
ac tions for six different flight altitudes. 

= TABLE 3 

- Geometrical Parameters for Observed Sea Area 

: 
r) 

2. 


e For the sea water we know that the average temperature 
in the geographical area of our interest is 13°C=286.15K for 
Midlatitude Summer and 18°C=291.15Kfor Tropical atmosphere 
[Ref. 6]. 
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For the sea water emissivity we chose the value of 0.96 
[Ref. 10:p.07]. 
For the calculation of the normal component of the sea 


radiance (W/cm?*Sr) in the range %=4.435ym to \.=5.00um we 
will use again the formulas from section A of chapter IV. 
This value must be multiplied by the water emissivity 
=0.96. For convenience we consider the observed area of sea 
as a point source. 

The normal component of the radiant intensity of this 
point source it will be: 
Louw Lye AREA 

And in the direction leading to the missile will be: 


Iw =Ton® cosa 
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V. DETECTOR MODELING AND CALCULATION OF S/N RATIO 


A. DETECTOR SYNTHESIS 


We already studied the geometry of the detection problem 
which shows us what the detector is required to do, and of 
course from this side we are interested in the most real- 
istic performance of the detector. Since we do not have 
available experimental data that could show us the perform- 
ance of the detector element we have to choose this 
according to the literature. 

From Figure 5.1 [Ref. 10] the possible choices for the 
detector include photovoltaic or photocontactive Indium 
Antimonide and Lead Selenide. The detailed data sheets for 
these detectors [Ref. 12:pp.365-368] show the relative 
advantages and disadvantages of each. When they are compared 
on the basis of D*, photovoltaic Indium Antimonide is 
clearly the best detector. A possible difficulty with this 
conclusion arises when the responsivities are compared and 
it is noted that value of Lead Selenide is 10 times higher 
than it is for the other two types of detector. In an ideal 
operating environment, that is, in which there are no extra- 
neous sources of noise, the designer is interested in the 
Signal - to - noise ratio of the detector, as described by 
D* and not the signal level, as described by responsivity. 
In advance we know that during the experimental trials 
photovoltaic Indium Antimonide had the better performance. 
Therefore we are led to the decision to use one element of 
photovoltaic Indium Antimonide, and we can see the detailed 
data for this material in Figure 5.2 [Ref. 12:p.368]. 

This one element will be placed -in the middle of an 


assumed two mirror system as in Figure 5.3, with 
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Figure 5.1 The D* of the Detectors in Each atte ber 
Window at the Most Common Operation Temperature OJ 
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Indium Antimonide (InSb) 


° 
rey 


Type: Photovoltaic. 

Element From.) x 0.} 1010 « 10 mm. 

sizes Rectangular elements from 

available: 0.1 mm wide. Circular elements 
with diameters from 0.1 to 10 mm. 

Time Less than | usec. 

constant: 

Dynamic 2-8 x 10* ohm. 

impedance: 

Responsivity: 2x 10*VW-'- 


D* (peak) 
D* (S00°K) 


D* (d, 1000) cm(Hz)* w-! 
ry 


~ 
roy 
3 


3 5 
Wavelength (microns) 


19 '2 


D* tpeak) cm(Hz) wo! 


180 
Temperature (*K) 


D*, 300°K 
background 


Noise voltage, ViHz)~ ‘ 


Element 1x1 mm | 
Bandwidth 1 Hz 


100 1000 


Frequency (Hz) Background photon flux (0-5 4y). 


photon cm~?sec ~! 


Figure 5.2 Data Sheet for Photovoltaic Indium 
Antimonide Detector [Ref. 13] 
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reflectivity 0.9 each mirror and total of two ba 
=0.9°0.9=0.81. The detector will be considered to be used 
with a two-mirror system designed so that the detector 
element will subtend area equal to the I.F.V. at a distance 
of PATH. Thus we have to determine the size of the sensitive 
element (detector). 


SEAF.CE 
TIELD 


SENOLDIVE 
CLEMENT 


Figure 5.3 System of Two Mirrors in the Receiver 


We know that the area of the detector A is given by: 
Ag=wef? 
where 
w=Instantaneous Field of View in (Sr) 


w= (SIDE*SIDE)/ (PATH) *=(2200¢2200)/(PATH)? (Sr) 
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f=equivelant focal length 
and for PATH we have to use the values from Table 2 which is 
function of the flight altitude. 

For the our purpose of work a very good and accurate 
approach to calculate the equivalent focal length of the 
system of two mirrors is the following. 

First we have to choose the diameter of the entrance 
aperture D,=l2cm. This dimension is good enough to have an 
aerodynamic shape in the nose of the missile. Thus from the 
Figure 5.3 we have: 

(AB)?=(AO)?+(0OB)?=72 and AB=8.4cm 
(AO)?=(AC)?+(CO)? or 
(CO)?=(A0)?-(AC)?=(AO)?-(AB/2)?= 
=(CO)?=6?-(8.49/2)?=17.98cm? 
and CO=4.24cm 
We can define the distance CO as the equivalent focal length 
£=4.24cm 
And 
Aa=wef?=(2200)7°(4.24)?/ (PATH)? 

Substituting the values of PATH for each flight altitude 
gives A and the square root of this gives us the side of 
the detector element.Table 4 demonstrates the dimensions of 
the detector element shows that with reasonable detector 
size the I.F.V. at about 1125 meters slant path range will 
include the hot region of the ship regardless of the flight 
altitude. 


B. SIGNAL VOLTAGE 


The general formula that gives the signal voltage is: 
+ 
Vs = [AsneR/ (PATH)? ] «Lf Tytaundl 
mt 
Where 
Ae :Area of entrance aperture 
Ao = WD? /4=113.1cm? 


ro :Optics reflectance 
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TABLE 4 


Dimensions of the Detector Element 


R :responsivity=2°10° (V/W) 
I, :Spectral radiant intensity of the source 
tT :Spectral atmospheric transmittance 
The above relation can be further simplified using 
a instead of the spectral transmittance an average one as we 
‘@) ; determined in Section C of Chapter III i.e. 0.4711 for 
. Midlatitude summer and 0.4259 for Tropical atmosphere. 
Taking Tala) outside of the integral then the quantity 
[tsa becomes simply the inband radiant intensity of the 
® source which has been calculated in section D of chapter 
III. Hence the signal voltage formula becomes: 
Vs=Ae re Rete Ly/ (PATH)? 2 
In order to determine the signal voltape due to all the 
e available point sources inside the I.F.V. we use the super- 
position principle: 
Vs=[AgreReTa/(PATH)?]e(IyslacIsslyele tie ely] 
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C. NOISE VOLTAGE 


The noise voltage is due to two main factors, namely the 


e background noise that is created from the background radiant 


intensity and the noise from the optical system. 


In Section E of Chapter III we determined the radiant 


intensity of the sea region that the detector observes and 


” we named Iw (W/Sr). The noise voltage of this reason will 
be: 


Vag? AS HP R OT" Lw/ (MK)? 


where MK is determined in Figure 4.3 and indicates 


€ the average path between the observed sea surface 


and the detector. 


optical system is given by: 


The noise signal due to the 
{Ref. 12]. 
. Vao=Re (Ag Af) # /D* 


where 


D¥:Detectivity with average value 9*10'? f[cem(Hz)?7?/wW]) in 


the wavelength range %=4.435pm to Ae5.00um 
Ag:Area of the detector as determined in previous section of 


the present chapter. 
R :Responsivity 


Qf :equivalent noise bandwidth 


C) The equivalent noise bandwidth is given by 
Af=1,2Te 


where 


Ta:dwell time (time required for the image of the target to 


e pass across the detector). 
and Tuewee 


where 


c:number of detector elements (1 for our caSe) 


e Q:Time rate of search. 


and 
Q-Q/Te 


where 


ae eS 
~ 


a] 


Q:size of search field in (Sr) 
Ge:frame time (time required to scan the entire search 
field). For our case we chose 1l0Omsec. 
Therefore . 
Af=1/2G=Q/2ewec=0/2ewecele 


PATHi 


Figure 5.4 Total Detector Search Field 


In Figure 5.4 you will find the geometrical representa- 
tion of the total search field which will have total length 
of 1000m and is scanned by the detector in =10msec. 

The total solid angle that the search field subtends at 
the detector is: 

Q= dQ 
dO=dx*SIDE/@’? 
£?=(PATH)?-x? 
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dQ=(SIDE/ (PATH)?-x?]edx 
QO=(SIDE) dx/ (PATH) ?-x? 


Since (PATH)? >500? always, from the tables we get: 


Q=(SIDE)/2(PATH) [1n((PATH+x) / (PATH-x))] 
= (SIDE) /2(PATH) [1n((PATH+500)/(PATH-500)) 

, - In( (PATH-500)/ (PATH+500) )] 

=((SIDE)/2(PATH) )*1n[(PATH+500) / (PATH-500) ]? 


Q=((SIDE)/ (PATH) )*1n[{(PATH+500)/ (PATH-500) J 


Therefore 


Af=(1/2(SIDE) *c*ts) *1n[ (PATH*+ 500) / (PATH-506)] 


The total noise voltage is 


Vn= Vug*t Vito 


And the signal to noise ratio is: 


Vs /Vn=Vs/ (Viet Veo) 


In figures 5.5 to 5.22 we can see the summary of the 
results of the above calculation for all the possible combi- 


nations of the relative parameters. 


va In Appendix D we can see the influence of sky condition 


to S/N ratio for the two extremes flight altitude i.e. 50 
(m) - 100 (m) and for all the combinations between atmos- 


pheric profile and paint emissivity. 
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VI. RESULTS AND DISCUSSION 


A. SHIP BODY TEMPERATURE 


‘Tables 5 to 18 of Appendix A show the effects of the 
ship internal temperature, atmospheric profile, ship body 
paint emissivity, sun elevation angle and sky condition on 
the ship body temperature. 

For. both horizontal and vertical faces higher paint 
emissivity and ship internal temperature result in higher 
ship body temperature. Higher internal temperature corre- 
sponds to less ship body heat losses towards the internal 
compartments. Moreover higher paint emissivity (absorp- 
tivity for gray paint) corresponds in higher solar energy 
absorption and, for the 300K ship body temperature level, 
the dominant mechanism of heat transfer is that of convec- 
tion. Hence the increase of the paint emissivity does not 
correspond to a significant increase of the energy radiated 
by the ship body, and the net result is an increase of the 
ship body temperature. 

The sky condition appears to be significant as well. The 
more clear it is the more. solar energy there is incident on 
the ship, resulting in an increase of its body temperature. 

As far as the atmospheric profile is concerned, the 
Tropical Atmosphere results in higher body temperature than 
that of Midlatitude Summer since the ambient temperature is 
higher and the ship body heat losses towards the environment 
are less. 

The sun elevation angle above the horizon affects 
vertical and horizontal faces in different ways. We know 
that the incident solar energy at sea level is an increasing 
function of the elevation angle. For an horizontal face the 
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temperature appears to be an increasing function of the 


elevation angle, the sine of which is used for the calcula- 
tion of the incident angle effect. For vertical faces the 
cosine of the elevation angle is used and the temperature 
appears a maximum in the vicinity of 50°. The explanation is 
that, up to this value, the increasing effect of the solar 
energy dominates the decreasing one of the cosine function. 
After the value of about 50° the first derivative of the 
solar energy function becomes small and the cosine behavior 
dominates so that the net result is a decreasing temperature 
function from 50° to 90°. 

AS we see in the ship tracked area geometric analysis is 
constituted from both vertical and horizontal faces, but the 
vertical ones have much larger projected surface areas. 
Therefore the radiant energy from the tracked area is due 
mainly to that of the vertical faces and as we will see from 
the S/N ratio results they appear maximum at a value of 
elevation angle close to 50°. 


B. DETECTOR S/N RATIO 


Figures 5.5 to 5.22 show the effects of the sky condi- 
tion, sun elevation angle and missile flight altitude on the 
detector S/N ratio. . 

In general we observe that the variation of the S/N 
ratio is similar to that of the ship body temperature. That 
means that the radiant energy from the ship body due to its 
temperature is the main thermal source for the detector. At 
this point, it is useful to discuss the influence of the 
Stack exit plane thermal radiation. From our calculation we 
know that the radiant intensity of the CO in the stack exit 
plane has a magnitude of the order of 10* (W/cm*Sr) and that 
of the ship body of the order of 10? to 10?. Moreover the 
stack exit plane area is much smaller than the total tracked 
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area and the stack exit radiant flux is significantly 
smaller than this of the ship body tracked surface. 

Sky condition and the paint emissivity affect S/N ratio 
similarly to the ship body temperature. The more clear is 
the sky and the higher is the paint emissivity the higher is 
the S/N ratio. 

The atmospheric profile affects S/N in a different way. 
While the Tropical Atmosphere gives higher ship body temper- 


ature, it results in lower S/N ratios comparing with the 
Midlatitude Summer. The explanation for that is that the 
increase of the ship body radiance, due to its temperature 


increase, is less than the background noise increase (higher 
sea temperature). Moreover the atmospheric transmittance in 
the 5um window for the Tropical Atmosphere is about 10% less 
than this for Midlatitude Summer. 

The missile flight altitude appears to be important as a 
parameter. The higher it is the lower is the S/N ratio and 
this decrease is higher between 50 to 70 (m) than between 70 
to 100 (m). The explanation is that although higher missile 
flight altitude corresponds to lower background area (sea), 
the path between missile head and ship becomes’ larger and 
its Square appears in the denominator of the signal voltage 
formula. Therefore the net result is to give lower values of 
S/N ratio. ; 

The influence of the sun elevation angle above the 
horizon is similar to that on the ship body temperature and 


it has been discussed earlier. 


C. RECOMMENDATIONS 


Based on the relative influence of each parameter which 


has been varied in our analysis we recommend: 


From the ship design point of view, lower internal 
temperature or, when it is not feasible, better wall 


insulation, especially in the area where significant heat 
generation exist. In addition to that, paint of low emis- 


Sivity and absortivity is desirable. 


From the missile point of view we recommend lower flight 
altitude. At this point we must notice that, lowering the 
flight altitude, we will reach a point where the background 


area (sea) will be very large, resulting in higher back- 


ground noise and the coast or even the sun itself could 
constitute background noise for the detector. Therefore the 
flight altitude must be determined, in a certain design, 
accounting for all these conflicting parameters. Another 
solution could be the use of higher flight altitude, that 
assures us low background area, and to improve the detector 
optics design that will balance the loss of the’ signal 
voltage due to the greater path length. 
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APPENDIX B 
TABLES FOR DETERMINATION OF I.F.V. 
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Path in Meters for Flight Altitude 
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TABLE 15 
Side of Instantaneous Field of View in Meters 
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TABLE 16 
Path in Meters for Flight Altitude 
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Side of Instantaneous Field of View in Meters 
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TABLE 18 
Path in Meters for Flight Altitude 
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TABLE 19 
Side of Instantaneous Field of View in Meters 
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TABLE 20 
Path in Meters for Flight Altitude= 80 m 
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TABLE 21 
Side of Instantaneous Field of View in Meters 
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TABLE 22 
Path in Meters for Flight Altitude 
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Side of Instantaneous Field of View in Meters 
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TABLE 25 
Data Used for Month October 


CARD 1 


Midlatitude Summer ; 

Slant Path Between Two Altitudes 
Prosran Execution in Transmittance Mode 
Midlatitude Summer T and P Profile - 
Midlatitude Summer Water Vapor, Profile 
Midlatitude Summer Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 


OOOOMMNONS!Y 


CARD 2 


3 Navy Maritime Extinction 

2 Fall - Winter . 

Q : No Stratospheric pacceroucd 

5 : Medium Continental Influence 

Q : No Cirrus 

Q Not Used 

0  : Meteorological Range(Km 

9: Current Wind Speed (m/s 

.09: 24 hours Average Wind Speed (m/s) 
0 : Rain Rate (mm/h) 


From Section G of Chapter II : Initial Altitude 
From Table 22 : Final Altitude 
-- :_ Not Necessary to Define 
From Table 22 : Path Length 
-- : Not Necessary to Define 
-- : Not Necessary to Define 
QO : Normal Operation 


2255 : Initial Frequency 
2000 : Final Frequency 
10 : Frequency Increment 
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TABLE 27 
Data Used for Month November 


CARD 1 

2 Midlatitude Summer ; 

2 Slant Path Between Two Altitudes 

Q Program Execution in Transmittance Mode 
2 Midlatitude Summer T and P Profile | 
2 : Midlatitude Summer Water Vapor Profile 
2 : Midlatitude Summer Ozone Profile 

0 Normal Operation of Program 

0 Normal Operation of Program 

0) Normal Operation of Program 

0 Normal Operation of Program 

CARD 2 

3 Bey Maritime Extinction 

2 Fali - Winter | 

Q : No Stratospheric Background 

5 : Medium Continental Influence 

Q : No Cirrus 

Q Not Used. 

10 Meteorological Range (Km) 

2.09: Current Wind Speed{m/ 8) 
2.09: 24 hours Ayer ees Wind Speed (m/s) 

QO : Rain Rate (mm/h) 


From Section G of Chapetr II : Initial Altitude 

Form Table 22 : Final Altitude 

-- : Not Necessary to Define 

From Table 22 : Path beneen 

-- : Not Necessary to Define 

-- : Not Necessary to Define 
Normal Operation 


2255 : Initial Frequency 
2000 : Final Frequency 
10 : Frequency Increment 


E19 


YIUOW ayy JO aseuaay [PIO] 


hOLh*O 


Spnitaiyd 


Z0Lh°0 2CLn'O EOLH*O EOLH'O 90Lh'D BDLn'O Gad Sh eaSAU 


E0Lh 0 EOLH"O nOLh 0 hOLH’O 


80zL4h°O TTL 0 (WA) 22079 


Gud) 02079 


(wi) RT0"O 


Z0Ln°O Z0Lh°O COLAO COLn'*O (U4) 3LO’O 
[ omy 
Cu 
LOLN‘O TOLh'O coe LOLNO CIN) 7 TAt9 n 
869h°0 369h'°0 eer 669h°O Z0Lh°O SOLh'O (UY) ONTO 
AQGNLIL IV : 
(WY ) (Wy) (Uy ) (UNt) (uy) (uN) EE LENT 
oot'o 060°9 Gane 0L0°0 030° 05079 2 
’ SOGnLILIV Re 
IWNIG : 


AJaqusAoN YAUOW TOF saouezqtTwsuer] advraay 


8c ATAVL 


TABLE 29 
Data Used for Month December 


Os = 


CARD 1 
MODEL 2 Midlatitude Summer ; 
ITYPE Z Slant Path Between Two Altitudes 
IEMSCT 0 Program Execution in Transmittance Mode 
M1 2 Midlatitude Summer T and P Profile | 
‘ M2 2 Midlatitude Summer Water Vapor Profile 
a M3 2 Midlatitude Summer Ozone Profile 
: IM 0 Normal Operation of Program 
NOPRT Q Normal Operation of Program 
TBOUND Q Normal Operation of Program 
: SALB 0 Normal Operation of Program 
CARD 2 
IHAZE 3 Navy Maritime Extinction 
ISEASN 2 Fall - Winter . 
IVULCN Q : No Stratospheric Background 
ICSTL 5 : Medium Continental Influence 
ICIR Q : No Cirrus 
. IVSA 0 Not Used. 
@ VIS 10 Meteorological Range 
- WSS 2.05: Current Wind Speed 
WHH 2.05: 24 hours Average Wind Speed (m/s) 
RAINRT O : Rain Rate (mm/h) 
CARD 3 
{ Hl : From Section G of Chapter II : Initial Altitude 
: From Table 22 : Final Altitude 
ANGLE : -- : Not Necessary to Define 
RANGE : From Table 22 : Path Length 
BETA : -- : Not Necessary to Define 
> -- : Not Necessary to Define 
LEN : O : Normal Operation 
OQ CARD 4 
ee Vil 2255 : Initial Frequency 
el V2 2000 : Final Frequency 
tae DV 10 : Frequency Increment 
& 


12) 


hOLh’O = UZUOW JUL Jo aBeuaAY [PIO] 


: ; oP yy 
Z0Lh*0 20Lhn°0 €0Lh*0 €0LH'0 902° 0 Pere ee 
hOLh*O hOoLnO SOLh*O SOLn‘O eee | (WA) 22070 
hOLh“O hOLH*O SOLH*O GOLH*O (UH) 


EOLH'O hOLH-O 


(Uy) 


oe 


Z0LH'O EOL’ O Cut) 


TOLh'O 20Lh' 0 (Unt) 


+ 
+ 


669°) OOLH*O IDEN OD (Uy) L200 


869t°O ( me (HN) HNO DO 


SONLILIV 
(uy) (uy ) (Yt) (ud) (Ut) (urd) ae 


0OT*O 060°0 080°90 GLO: |: 090°0 0S0°0 
SqNLILV 


‘IVNIG 
ATaquaceq YyAUOW AOF Saoue AtTwSUeI] BdeIZAY 


O€ ATAVL 


TABLE 31 
Data Used for Month January 


CARD 1 


Midlatitude Summer ; 

Slant Path Between Two Altitudes 
program Execution in Transmittance Mode 
Midlatitude Summer T and P Profile | 
Midlatitude Summer Water Vapor Profile 
Midlatitude Summer Ozone Profile 

Normal Qperation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 


ODOOOMNNOPNY 


CARD 2 


3 Navy Maritime Extinction 
2 Fall - Winter 
QO : No er peegsphertc Background 
5 : Medium Continental influence 
Q : No Cirrus 
0 Not Used 
15_(: Meteorological Range 
2.06: Current Wind Speed_ 
2.06: 24 hours Average Wind Speed (m/s) 
QO : Rain Rate (mm/h) 


From Section G of Chapter II : Initial Altitude 
From Table 22 : Final Altitude 

-- :_Not Necessary to Define 

From Table 22 : Path Length 

-- +: Not Necessary to Define 

-- : Not Necessary to Define 


O : Normal Operation 
CARD 4 
2255 : Initial Frequency 
2000 : Final Frequency 

10 : Frequency Increment 
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TABLE 33 
Data Used for Month February 


CARD 1 


Midlatiitude Summer ; 
Slant Path Between Two Altitudes 
: Program Execution in Transmittance Mode 
: Midlatitude Summer T and P Profile | 
: Midlatitude Summer Water Vapor Profile 
Midlatitude Summer Ozone Profile 
Normal Operation of Program 
Normal Operation of Program 
Normal Operation of Program 
Normal Operation of Program 


DOOCOMMNO!S)> 


CARD 2 


Navy Maritime Extinction 
Fall - Winter . 
No Stratospheric Background 
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Meteorological Range (Km) 

22: Current Wind Speed (m/s 

.22: 24 hours BY ET Oe Wind Speed (m/s) 
O : Rain Rate (mm/h) 


CARD 3 


Ner 


From Section G of Chapter II : Initial Altitude 
From Table 22 ; Final Altitude 

-- : Not Necessary to Define 

From Table 22 : Path beveen 7 

-- : Not Necessary to Define 

-- : Not Necessary to Define 

0 : Normal Operation 


2255 : Initial Frequency 
2000 : Final Frequency 
10 : Frequency Increment 


125 


8O0Lh°0 


OTL’ 0 


OtLh'O 


60Lh°0 


B0Lh°0 


LOLh*O 


S0zLh°0 


hOLh’O 


0oT*O 


OTLh’O 


802h'0 


OTLh'O 


cu) 


080°0 


60Lh'O 


OTLh'O 
60Lh°O 


od 


Cay) 


020°0 


ann nn Oe 


URUOW Buy JO aBeuaaay 


80Lh°0 


(uy) 


090°0 


Azeniqay yyuoW tox savouej.RTusuesy 


7€ DVL 


Teo] 


STLh*O 
LTLhO 


pare | 
STL O 
J sane | 


Tteh'O 


(ut) 


aZde1oAy 


ePpnartily 
Jad ageuaay 


(uy) 220°0 
(Wy) T20°O 


(44) B8T0°O 


(YX) 3T0°0O 


(wy) ZTO°9 


(YN) Lo00'd 


(44)  =f00°9O 


SGNLILTV 
TWILINI 


AGNLILTV 
TWNId 


TABLE 35 
Data Used for Month March 


CARD 1 


Midlatitude Summer ; 

Slant Path Between Two Altitudes 
Program Execution in Transmittance Mode 
Midlatitude Summer T and P Profile | 
Midlatitude Summer Water Vapor Profile 
Midlatitude Summer Ozone Profile 

Normal Qperation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 


OCOOOMNMMON)N 


3 Navy Maritime Extinction 

2 Fall - Winter . 

0 : No perergsr Jette Background 
5 : Medium Continental Influence 
p : No Cirrus 

5 


1 Meteorological Range (Km) 

2.13: Current Wind Speed (m/s) 

2.13: 24 hours Average Wind Speed (m/s) 
QO : Rain Rate (mm/h) 

CARD 3 


From Section G of Chapter II : Initial Altitude 
From Table 22 : Final Altitude © 

-- :_Not Necessary to Define 

From Table 22 : Path Length 

-- +: Not Necessary to Define 

-- : Not Necessary to Define 


QO : Normal Operation 
CARD 4 
2255 : Initial Frequency 
2000 : Final Frequency 

10 : Frequency Increment 
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TABLE 37 
Data Used for Month April 


CARD 1 
MODEL 2 : Midlatitude Summer . 
ITYP 2 : Slant Path Between Two Altitudes 
LEMSCT Q : Program Execution in Transmittance Mode 
M1 2 : Midlatitude Summer T and P Profile | 
M2 2 : Midlatitude Summer Water Vapor Profile 
M3 2 : Midlatitude Summer Ozone Profile 
IM 0 Normal Operation of Program 
NOPRT 0 Normal Operation of Program 
TBOUND 0 Normal Operation of Program 
ALB 0 Normal Operation of Program 
CARD 2 
IHAZE ae Navy Maritime Extinction 
ISEASN 2 : Fall - Winter | 
IVULCN Q : No SET ot osPuerse Background 
ICSTL 5 : Medium Continental Influence 
ICIR Q : No Cirrus 
IVSA Q : Not Used 
VIS 15__: Meteorological Range (Km) 
WSS : 2.12: Current Wind Speed (m/s) 
WHH : 2.12: 24 hours Aver ane Wind Speed (m/s) 
RAINRT : O : Rain Rate (mm/h) 
CARD 3 
H1 : From Section G of Chapter II : Initial Altitude 
H2 : From Table 22 : Final Altitude 
ANGLE : -- ?_Not Necessary to Define 
RANGE : From Table 22 : Path Length 
BETA : --  : Not Necessary to Define 
: -- +: Not Necessary to Define 
LEN : QO : Normal Operation 
CARD 4 
Vi 2255 : Imitial Frequency 
V2 2000 : Final Frequency 
DV 10 : Frequency Increment 
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TABLE 39 
Data Used for Month May 


CARD 1 
1 Tropical Atmosphere ; 
2 Slant Path Between Two Altitudes 
0 Program_Execution in Transmittance Mode 
1 Tropical Temperature and Presure Profile 
1 Tropical Water Vapor Profile 
1 Tropical Ozone Profile 
0 Normal Qperation of Program 
Q Normal Operation of Program 
0 Normal Operation of Program 
Q Normal Operation of Program 
CARD 2 
3 Navy Maritime Extinction 
1 Spring - Summer 
0 No era cospheric Background 
5 Medium Continental Influence 
0 No Cirrus 
0 Not Used. 
15 Meteorological Range pea) 
2.02: Current Wind Speed (m/s) 
2.02: 24 hours Averaee Wind Speed (m/s) 
0 Rain Rate (mm/h) 
CARD 3 
From Section G of Chapter II : Initial Altitude 
From Table 22 : Final Altitude 


CARD 4 


0 


Not Necessary to Define 


From Table 22 :;: Path Length 


Not Necessary to Define 
Not Necessary to Define 
Normal Operation 


Initial Frequency 
Final Frequency 
Frequency Increment 
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TABLE 41 
Data Used for Month June 


CARD 1 
MODEL 1 Tropical Atmosphere ; 
ITYPE 2 Slant Path Between Two Altitudes 
IEMSCT 0 Program Execution in Transmittance Mode 
M1 1 Tropical Temperature and Presure Profile 
\ M2 1 Tropical Water Vapor Profile 
M3 1 Tropical Ozone Profile 
IM 0 Normal Operation of Program 
NOPRT 0 Normal Operation of Program 
TBOUND Q : Normal Qperation of Program 
SALB Q : Normal Operation of Program 
j CARD 2 
i 
IHAZE 3 Navy Maritime Extinction 
ISEASN 1 Spring - Summer 
IVULCN 0 No Stratospheric peeneroune 
ICSTL 5 Medium Continental Influence 
ICIR Q No Cirrus 
: IVSA 0 Not Used 
@ VIS 12 Meteorological Range (Km) 
WSS 2.15: Current Wind Speed (m/s) 
WHH 12,15: 24 hours et a Wind Speed (m/s) 
RAINRT : O : Rain Rate (mm/h) 
CARD 3 
{ H1 : From Section G of Chapter II : Initial Altitude 
H2 : From Table 22 : Final Altitude 
ANGLE : --  : Not Necessary to Define 
RANGE : From Table 22 : Path Length 
BETA > -- +? Not Necessary to Define 
RO totes, 93, NOE Dp aceon Abe Define 
LEN : QO : Normal Operation 
vr) CARD 4 
: V1 : 2255 : Initial Frequency 
V2 : 2000 : Final Frequency 
DV : 10 : Frequency Increment 
© 
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Pale TABLE 43 


Pi Data Used for Month July 
me CARD 1 
MODEL 1 Tropical Atmosphere ; 
ITYPE 2 Slant Path Between Two Altitudes 
LEMSCT 0 Program Execution in Transmittance Mode 
M1 1 Tropical Temperature and Presure Profile 
M2 1 Tropical Water Vapor Profile 
M3 1 Tropical Ozone Profile 
M 0 Normal Operation of Program 
OPRT 0 Norma] Operation of Program 
TBOUND Q Normal Operation of Program 
ALB 0 Normal Operation of Program 


P | “CARD 2 


IHAZE 3 Navy Maritime Extinction 
ISEASN 1 Spring - Summer 
IVULCN Q : No et hateseuerte Background 
ICSTL 5 : Medium Continental Influence 
ICIR Q : No Cirrus 

. IVSA Q Not Used 

e VIS 12. _: Meteorological Range en) 
WSS 2.15: Current Wind Speed (m/s) 
WHH 2.15: 24 hours AveTaee Wind Speed (m/s) 
RAINRT 0 Rain Rate (mm/h) 

CARD 3 

{ Hl : From Section G of Chapter II : Initial Altitude 

ra : From Table 22 : Final Altitude 
ANGLE : -- : Not Necessary to Define 
RANGE : From Table 22 : Path Length 
BETA : -- : Not Necessary to Define 

>: -: : Not Necessary to Define 

LEN : QO : Normal Operation 

O CARD 4 
V1 : 2255 : Initial Frequency 
V2 >: 2000 : Final Frequency 
DV : 10 : Frequency Increment 
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TABLE 45 
Data Used for Month August 


CARD 1 
MODEL 1 Tropical Atmosphere : 
ITYPE 2 Slant Path Between Two Altitudes 
LEMSCT 0 Program Execution in Transmittance Mode 
M1 5 Tropical Temperature and Presure Profile 
M2 1 Tropical Water Vapor Profile 
M3 1 Tropical Ozone Profile 
9) Normal Operation of Program 
PRT Q Normal Operation of Program 
TBOUND Q Normal Operation of Program 
ALB 0 Normal Operation of Program 
CARD 2 
IHAZE 3 Navy Maritime Extinction 
ISEASN 1 Spring - Summer 
IVULCN QO : No Stratospheric Background 
ICSTL 5 : Medium Continental Influence 
ICIR Q : No Cirrus 
IVSA 0 Not Used 
VIS 10 Meteorological Range (Km) 
WSS : 2.20: Current Wind Speed (m/s) 
WHH : 2.20: 24 hours Average Wind Speed (m/s) 
RAINRT : QO : Rain Rate (mm/h) 
CARD 3 
Hl : From Section G of Chapter II ; Initial Altitude 
: From Table 22 : Final Altitude 
ANGLE -- : Not Necessary to Define 
NG From Table 22 : Path Length 
BETA > -- : Not Necessary to Define 
RO > -> : Not Necessary to Define 
LEN >: QO : Normal Operation 
CARD 4 
V1 2255 : Initial Frequency 
V2 2000 : Final Frequency 
DV 10 : Frequency Increment 
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TABLE 47 
Data Used for Month September 


CARD 1 


Tropical Atmosphere ; 

Slant Path Between Two Altitudes 

Program Execution in Transmittance Mode 
Tropical Temperature and Presure Profile 
Tropical Water Vapor Profile 

Tropical Ozone Profile 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 

Normal Operation of Program 


OOOCrrFrFOhM sr 


Navy Maritime Extinction 

Spring - Summer 

No Stratospheric Background 
Medium Continental Influence 

No Cirrus 

Not Used 

Meteorological Range (Km) 

Current Wind Speed (m/s) 

24 hours Average Wind Speed (m/s) 
Rain Rate (mm/h) 
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From Section G of Chapter IL : Initial Altitude 
From Table 22 : Final Altitude 
~- i: Not Necessary to. Define 
From Table 22 : Path Length 
~- : Not Necessary to Define 
~>- : Not Necessary to Define 
Normal Operation 


Initial Frequency 
Final Frequency 
Frequency Increment 
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APPENDIX D 
PLOTS OF S/N RATIO 
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APPENDIX E 
COMPUTER PROGRAMS 
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C77= (JEIGAT= #37) 4eaLa 
aL? (3 L)=KI7NSC 
41 CONTINGE 
42 CONTINUE 
pes 
¢ CALCULATION OF SHIP SIGNAL VOLIac2 
DO 50 _K=1,3 
oes te bate, L) +aL2 (KL) eXL3 (XL) OSTA (K 315 jK I 
= K +ale Ly ed 4 ¢ pt L) ¢ be 
ie RIT IR Lh) TRA US a RE EC wor Tek DORE (EAT BEEK TN)” 
49 CONTINIE 
é 50 CONLING: 
€ CALCULATICN UF NCISE VOLiaGE 
usldst=sLoeestge/ (fATH@vATAD 
AD=JSINST#ECCL EGU 
QHTOT> (S10E/ eA chy ALU | (ca Zhe xtOT) / (2ATH-UTCT) | 
DERE G=UATOT (2. SUA LNST*EN®TLACE) 
VuoPT-2eSECus{(adeGESE 4) ** 23) /LELECT . 
VB2ddL a ®AREASFLRAAS HRESCCH SILA DEAE ESTY (DQATHS*PATHS) 
: UNTOL=VNCPL eva 
€ CALCULATION UF SIGNAL TC NOISE waite 
30 60 421,45 : 
00 39 L2i,91 
S¥P (K LF =20. *ALOGII (VS (KL) /VATIT) 
$9 CONTINUE 
60 CONTINUE 
E 
€ TASULATION OF Thd DETECTOR SIGhAL Tu NCISE 3aTIC 
aRITE (0,70) 
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Steere eft sat HEIGHL=50 MbigeS aND SILE=Z1.9d YEGERS"//) 
nN FORMAT (Sy!) att. acca! 3h, *OUGHSCIEED SUN',3%,'SUN «IT LIGHL CLOO 
10°, 5k, (SON ali HEAVE SIGS CLOJDS'//) 
00°73 121,91 
WRITE (Dd, 12) I (Ld y (HS (KL) 21, 
12 SOBIR S| eer ee othe F1003,76,F10.5, 138,710. 3) 
73° CONTINU 
STce 
BND 
SENTRY 
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